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We experimentally examine the nonclassical character of a class of non-Gaussian states known as phase- 
diffused squeezed states. These states may show no squeezing effect at all, and therefore provide an interesting 
example to test nonclassicality criteria. The characteristic function of the Glauber-Sudarshan representation 
(P function) proves to be a powerful tool to detect nonclassicality. Using this criterion we find that phase- 
diffused squeezed states are always nonclassical, even if the squeezing effect vanishes. Testing other criteria of 
nonclassicality based on higher-order squeezing and the positive semidefinitness of special matrices of normally 
ordered moments, it is found that these criteria fail to reveal the nonclassicality for some of the prepared phase- 
diffused squeezed states. 

PACS numbers: 03.65.Wj, 42.50.Dv, 42.50.Xa 



I. INTRODUCTION 

The definition of nonclassicafity of a quantum state of the 
harmonic oscillator is closely connected to the coherent states. 
These are the eigenstates of the annihilation operator, a\a) = 
a I a) , where the complex number a defines the amplitude and 
phase of the field Glauber and Sudarshan showed that the 
density operator of an arbitrary optical quantum state can be 
formally written as a statistical mixture of coherent states. 
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where the Glauber-Sudarshan representation P{a) plays the 
role of the probability distribution of coherent states |0, Ht]. 
However, in quantum optics P{a) often violates the proper- 
ties of a probability density. Hence, a state is referred to as 
nonclassical if its P function does not exhibit the properties 
of a classical probability density 

Only recently, nonclassicality of experimentally generated 
states has been demonstrated by means of this definition jsj]. 
In many cases, however, the P functions of nonclassical states 
are highly singular, such that they cannot be reconstructed 
from the measured experimental data. This is the case for 
squeezed states, having a quadrature variance of less than the 
quadrature variance of the vacuum state. For instance, the P 
function of a squeezed vacuum state with quadrature variances 
Vx and Vp (we assume that Vx < 1, where unity represents the 
normalized vacuum noise) may be formally written as 







Sia). (2) 



This quantity cannot be understood as a well-behaved func- 
tion. 

Phase-diffused squeezed states define an interesting class 
of states with a, in general, not accessible P function. In 
very recent experiments these states were used to demonstrate 
purification and distillation for continuous variable quantum 



information protocols ^ . Phase diffused squeezed states 
are a mixture of squeezed (vacuum) states with a stochasti- 
cally distributed phase. They are related to a realistic decoher- 
ence process and may be produced from pure squeezed states 
in a phase noisy transmission channel. They reveal a non- 
Gaussian noise distribution, have a positive Wigner function, 
and, for strong phase noise, may show no squeezing effect at 
the level of second moments of the quadrature operators. 

In this paper, we use phase-diffused squeezed states in or- 
der to experimentally test nonclassicality criteria for the case 
where the P function cannot directly be reconstructed from 
the homodyne detector quadrature data. First, we concentrate 
on the characteristic function of the P function, which is al- 
ways well-behaved, and investigate the criterion proposed in 
Ht]. Second, we examine moments of the quadrature operator 
and search for higher-order squeezing . Third, we check a 
hierarchy of criteria based on iiormally ordered moments, as 
suggested in ifioll . We find that the characteristic function of 
the P function outperforms the other criteria of nonclassical- 
ity. 

Let us consider a statistical mixture of squeezed states, each 
described by a Wigner function [il lil 
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where = xcos{if) +psin((y9) and — ~xsm{Lp) + 
pcos(iy9) are the quadrature variables, rotated around an angle 
<f, and Vx: ^re the variances of both quadratures x^p^p^p, 
satisfying the Heisenberg uncertainty relation VxVp > 1. Let 
p{(p) denote the statistical distribution of the phase fluctua- 
tions, then the Wigner function of the mixed state reads 



W{x,p)^ p{ip)Wsv{x,p;ip)dip. 



(4) 



In our examination of nonclassicality, the characteristic func- 
tion $(/?) of the P function plays a decisive role. It is con- 
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Table I: Parameter of the examined states, 
nected to the Wigner function via Fourier transform 

For a squeezed state, as defined by Eq. (O, we find 

= g^[l-VxCos2(arg(/3)-<p)-ypSin^(arg(/3)-<p)]^ 

(6) 

The characteristic function for the mixed state is given in close 
analogy to Eq. (|4|, 



(7) 



In our experiment we generated phase-diffused squeezed 
vacuum states with varying strengths of the phase noise. The 
phase noise was chosen to be distributed according to a zero 
mean Gaussian and could therefore be completely character- 
ized by the standard deviation. A summary of states generated 
is given in Table HI The undisturbed squeezed vacuum states 
had quadrature variances Vx = 0.36 and Vp = 5.28. For the 
strongest phase noise we used a flat distribution with a width 
of 720°, which is labeled with cr = oo in Table U We also 
listed the minimum quadrature variance of each state. 
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to show that the states with a = 6.3° and a = 12.6° are still 
squeezed, but the squeezing vanishes at cr = 22.2°. Hence, 
one cannot decide about the nonclassicality of the last two 
states by examination of the quadrature variance. 



II. EXPERIMENTAL SET-UP 

The squeezed states were generated by a degenerate opti- 
cal parametric amplifier (OPA). The OPA consisted of a type- 
I non-critically phase matched second order nonlinear crys- 
tal (7% MgiLiNbOs) inside a standing wave optical resonator 
with a line width of 25 MHz. The OPA process was contin- 
uously pumped by 50 mW of second harmonic light yielding 
a classical power amplification factor of six. Both the length 
(resonance frequency) of the resonator as well as the orienta- 
tion of the squeezing ellipse were stably controlled by elec- 
tronic servo loops. With this setup we directly measured a 
minimal squeezed variance of -4.5 dB and an anti-squeezed 
variance of H-7.2 dB with respect to the unity vacuum variance. 
From these measurements we inferred an overall efficiency of 
75% and an initial squeezing factor of -8.2 dB. 

The squeezed field propagated in free space from the OPA 
passing high-reflection mirrors, two of which were moved by 
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Figure 1: Simplified sketch of the experimental setup. MC: spa- 
tial mode cleaner, LO: local oscillator, OPA: squeezed light source, 
BHD: balanced homodyne detector, PZT: piezo-electrically actuated 
mirror. PZTl was used to control the average phase and PZT2 ap- 
plied the phase noise. 



piezo-electric transducers (PZT). One (PZTl, Fig.[T]i was used 
to control the average phase of the squeezed field. The other 
(PZT2, Fig. [T]i was driven by a quasi-random voltage to ap- 
ply the phase diffusion. This voltage was generated by a high 
quality PC sound card connected to an appropriate amplifier. 
The sound card played back a previously generated sound file 
which was carefully designed to meet the desired shape of its 
frequency spectrum and its histogram. The former covered 
the flat part of the frequency response of the PZT but the fre- 
quency band of any control loop, the latter was chosen to be 
Gaussian for the partial phase diffusion and had to be abso- 
lutely flat in the totally randomized case. 

Balanced homodyne detection (BHD) was used to mea- 
sure the quadrature amplitude of the phase-diffused squeezed 
field. The visibility of the squeezed beam and the spatially 
filtered (MC, Fig.[T]) local oscillator was 98.9% and was lim- 
ited by OPA crystal inhomogeneities. The average quadrature 
phase of the BHD was servo loop controlled except for the 
total phase randomization where no mean phase exists. The 
signals of the two individual BHD-photodetectors were elec- 
tronically mixed down at 7 MHz and low pass filtered with a 
bandwidth of 400 kHz to address a modulation mode showing 
good squeezing and a high dark noise clearance of the order 
of 20 dB. The resulting signals were fed into a PC based data 
acquisition system and sampled with one million samples per 
second and 14 bit resolution. For a more detailed description 
of the main parts of the setup we refer to |0, 0] • 



III. NONCLASSICALITY IN TERMS OF THE 
CHARACTERISTIC FUNCTION 

A. Experimental demonstration 

First, let us consider a sufficient criterion proposed in Ht): A 
state is nonclassical if the characteristic function $(/3) of the 
P function exceeds the characteristic function of the vacuum 
at some point (3: 



3(3 : > I*, 



(9) 



Note that this condition represents the lowest order of a hi- 
erarchy of conditions which completely characterize the non- 
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Figure 2: Characteristic functions of different pfiase-diffused 
squeezed vacuum states. Tfie sfiaded area corresponds to one stan- 
dard deviation, it is added to tfie nonclassical boundary <E>vac = 1- 
Take note of tfie logaritiimic scale. 
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classicality lll2ll . The function can be obtained by |ll3| 



^{13) = (: b{p) :) = ^e^l/3|-(V2-arg(/3))\^g|/3|V2^ 



(10) 



where D{j3) is the displacement operator. Since we only con- 
sider a single quadrature, we may neglect the arguments of (3 
and a;. The expectation value on the right hand side of Eq. ( fTOl i 
represents the characteristic function of the quadrature. It can 
be estimated from the sample of N measured quadrature val- 



ues{a;j}™^^ via (cf. il4ll) 
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Inserting Eq. ( fTTT i into ( fTOl ). we obtain an estimation <l>(/3) of 
<I'(/3). The variance of this quantity can be estimated as 



{¥(/?)} = l[e>«l^-|¥(/3)|^ 



(12) 



For each state we have recorded = 10' data points. The 
resulting characteristic functions are shown in Fig.|2] We only 
concentrate on the quadrature where the variance of the state 
is minimum. The shaded area corresponds to the magnitude 
of one standard deviation, it is added to the nonclassicality 
border I $ vac I = 1- In order to demonstrate that a state sat- 
isfies the nonclassicality criterion, Eq. (|9]l, with a significance 
of s standard deviations, we have to check if the characteristic 
function satisfies the inequality 



W)\ > 1 



{$(/?)} 



(13) 



at least at one point /?. It is clearly seen that all recorded states 
satisfy this lowest-order condition for nonclassicality with a 
high significance of s > 10. Hence, we directly observe sig- 
natures of nonclassicality in the characteristic functions. 

We note that Fig. |2] reveals that the P functions of these 
states are highly singular. This is due to the fact that condi- 
tion (|9]l is satisfied for all (3 with large modulus, indicating that 



$(/?) is not integrable. Therefore, we cannot expect P{a) to 
be a well-behaved function. However, if does not satisfy 
Eq. (|9]l, one may be able to compute its Fourier transform and 
check nonclassicality of the state based on the failure of the P 
function to be nonnegative, see |5ij]. 



B. Theoretical generalization of the results 

Whereas in the experiment only states with Gaussian phase 
noise are investigated, we now prove that phase-diffused 
squeezed vacuum states always fulfill condition (|9]l, indiffer- 
ent of the phase distribution p{(p). Although this criterion is 
not necessarily fulfilled by an arbitrary nonclassical state, it 
turns out to be sufficient for any state with a characteristic 
function of the form (|7]i. 

First, we note that <l>sv(/3; f) is 7r-periodic in the angle (p. 
Without any loss of generality, we can assume that the func- 
tion p{(p) in Eq. (|7]i can be regarded as a probability density 
over the interval [0, tt], i.e. p{f) > and p{(fi)dip ~ 1. 
For technical reasons, we may further assume that p{(p) is tt- 
periodic as the characteristic function of squeezed vacuum, 
such that the domain of integration in Eq. ^ may be any in- 
terval of length TT. 

It is clear that the average value of the function p(</?) equals 
to TT^^ on any interval of length tt. As an intermediate result, 
we show that for any positive integer n there exists an interval 
[■^0 ^ 5^ ' '/'o + on which the average of p{(p) is greater or 
equal than tt^^: 



1 TT 1 

p{ip + (pa)dip > = -. (14) 

IT n n 



Let us assume that this claim is wrong. Then we can find a 
certain n such that 



p{ip + ipo)d(p < - 
n 



(15) 



holds true. With this assumption we can decompose the inter- 
val [0, tt] into a set of n successive intervals of the length ^ 
and calculate 

rp{^)d^ = E + 0- - m < E ^ - 1^ 

(16) 

which is in clear contradiction to the normalization condition 
of p{(p). Consequently, our assumption is wrong, and the de- 
sired statement is verified. 

Now let us look at the characteristic function defined in 
Eqs. (|6]l and (|7]i. We assume that the z-quadrature is the 
squeezed one, so Vx < 1 < Vp. Hence, we can find a positive 
integer n such that 



" 2n ' 2n -I 



1 - Vx 008^(^3) - Vp sin2(^) > (17) 



holds. Simultaneously, we can find a real number ipo such 
that the inequality ( fT4b is fulfilled, and choose /3 = \f3\e'-f°. 
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Making use of the simple relations 

$sv(/3; - *sv(/3e-*'^; 0), $,v(/3; 0) = $sv(/3*; 0), 

(18) 

and taking into consideration that the domain of integration in 
Eq. (|7]i may be an arbitrary interval of length tt, we calculate 



ripa + -2 

e'^-) = / ^ p{ip)^,,i\(3\e'^^';ip)dip 
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P(^)1>sv(|/3| 



,l(¥)o-¥j). 



;0) 



;0)d^. (19) 



Since the integrand is nonnegative, a diminution of the domain 
of integration decreases the value of the integral: 



p(¥. + ^o)$sv(|/3|e''";0)d<^. (20) 



Because n is chosen such that Eq. ( fTTb is fulfilled for each </? 
in the domain of integration, there exists an e > such that 



(21) 



Consequently, we have a lower bound for Eq. (|2Q] t: 

_ZL 

l^d/^le^'^^iO)! > e'\^\"/' p{^ + ^o)d^ (22) 

Additionally, we have chosen (/?o such that Eq. ( fT4l ) is fulfilled. 
Therefore, 



e|/3|V2 



(23) 



Obviously, the characteristic function <i>(/3) is not bounded, 
independently of the phase noise distribution p{ip). Conse- 
quently, we are always able to prove nonclassicality by means 
of Eq. (|9]), which is the most simple criterion among a neces- 
sary and sufficient hierarchy lfT2ll . 



IV. NONCLASSICALITY IN TERMS OF MOMENTS 

Besides the signatures of nonclassicality in terms of the 
characteristic function, different criteria for demonstrating 
nonclassicality are known. Since we measured time series of 
individual quadrature values we can calculate all higher-order 
moments of the quadrature operator x. 



A. Hong-Mandel squeezing 

We begin by examining higher-order squeezing as proposed 
by Hong and Mandel @] . To this end, we calculate the degree 



of 2n-th order squeezing 

((Ai)2») 
(2n- 1)!! 



92n 



1, 



n e N 



(24) 



where Aa; = x — (x). The moments can be estimated from 
the sample of quadrature data quite naturally by replacing ex- 
pectation values by their arithmetic means. It is sufficient to 
verify nonclassicality if at least one of the q2n is negative. 

Table shows the degree of squeezing for different or- 
ders and different phase noise strengths. Obviously, only if 
the lowest-order parameter q2 is negative, then the parame- 
ter of higher order can also be negative. Hence, we may 
only observe higher-order squeezing if the state already shows 
standard squeezing. The investigation of the degree of higher- 
order squeezing does not extend the range of detection of non- 
classicality of phase-diffused squeezed states. This is not sur- 
prising, since it can be shown that for Gaussian states (O the 
degree of squeezing is given by 



q27i{v) = [Vx cos{lp) + ypsiii(^)]" - 1. 



(25) 



For these states, squeezing always implies higher-order 
squeezing and vice versa i^. Phase diffusion can only smooth 
out the phase dependence of the moments and diminish the 
nonclassical effect. 



B. Matrices of normally ordered moments 

Normally ordered moments of the quadrature operator x 
can be estimated from measured data points {xj}jLi via ap- 
propriate sampling relations (see AppendixlAli. 



1 



N 



2fe/2j\r 



V2 



(26) 



where Hk{x) are the Hermite polynomials. With these mo- 
ments at hand, we can examine the nonclassicahty criterion 
of Agarwal fioll. It has been shown that a state is nonclassical 
if at least one of the matrices 



i-.x:) {-.x':) 



X :) 

2 



( 



X' -.)\ 



X 



(27) 



is not positive semidefinite. This can be verified by showing 
that at least one of the principal minors of such a matrix is 
negative ifTsll . However, to this end we had to check up to 
2' — 1 principal minors for each matrix AI^''\ which is a com- 
putationally expensive task. 

Here we use the fact that the existence of a negative eigen- 
value of M^') demonstrates the violation of positive semidef- 
initeness. Therefore, we determine 



.(/) _ . x^M^'^x 

^min - -JT^ 



(28) 
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0.0 
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-0.6362(1 ± 0.3%) 
-0.5717(1 ±0.04%) 
-0.4060(1 ±0.08%) 
0.0196(1 ±3.2%) 
1.908(1 ±0.09%) 


-0.8667(1 ±0.16%) 
-0.8090(1 ± 0.03%) 
-0.5509(1 ±0.15%) 
0.6864(1 ± 0.53%) 
10.68(1 ± 0.16%) 


-0.9506(1 ± 0.12%) 
-0.9102(1 ± 0.03%) 
-0.5459(1 ± 0.60%) 
2.982(1 ± 0.84%) 
51.72(1 ±0.32%) 


-0.9813(1 ±0.09%) 
-0.9549(1 ± 0.04%) 
-0.3852(1 ±4.2%) 
10.61(1 ± 1.7%) 
249.6(1 ± 0.65%) 


-0.9927(1 ±0.07%) 
-0.9754(1 ±0.06%) 
0.0798(1 ± 95%) 
37.27(1 ± 3.3%) 
1222(1 ± 1.23%) 



Table II: Degree of squeezing q2n for different orders 2n and standard deviations of phase noise a. 



via a conjugate gradient algorithm, see, e.g., fldn . It can be 
shown that X^^j^^^ equals to the minimum eigenvalue of il/'''. 
In this way, we only need to calculate one quantity per matrix 
to examine its definiteness. Its standard deviation is deter- 
mined by using a bootstrap method: We generate new quadra- 
ture data, distributed as the experimentally measured quadra- 
tures, 100 times to obtain a statistical sample of eigenvalues, 
which gives the standard deviation, cf . IitIi . 

The experimental results are shown in Table |III] We only 
consider matrices of even dimension, since we noted that the 
minimum eigenvalues of M^^") and A/(2"+i) are equal. This 
may be due to the fact that odd moments of squeezed vacuum 
states vanish, giving the matrices a special structure. We ob- 
serve that all matrices, which belong to the states showing 
squeezing, have significantly negative eigenvalues. Hence, 
their nonclassical character can be directly observed in the 
sign of the smallest eigenvalue. Furthermore, for the states 
with (T > 22.2° the matrix M*^^' is positive semidefinite, since 
this directly corresponds to the absence of quadrature squeez- 
ing. However, for the state with a = 22.2° the matrices A/*^'^ 
with I > 4 possess a negative eigenvalue. Therefore, Agar- 
wal's criterion extends the range of detection of nonclassical- 
ity. Only for the completely phase-diffused state, are we not 
able to prove nonclassicality by this method. For this state, 
the effect might appear in higher-dimensional matrices, but 
the statistical uncertainty might hide the effect. 



ments. 

Eventually, we note that the evaluation of the statistical sig- 
nificance of nonclassical effects is much easier in terms of 
characteristic functions, since Eq. (fT2l i provides a simple rela- 
tion between the variance and the value of Testing the 
definiteness of matrices of moments requires complex non- 
linear procedures, for instance the calculation of the smallest 
eigenvalue or the principal minors. This leads to complica- 
tions in the estimation of the statistical significance, which 
increase the computational effort. For matrices of high orders 
the resulting errors are large and they may hide the sought 
nonclassical effects. 
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Appendix A: SAMPLING FORMULA FOR NORMALLY 
ORDERED MOMENTS 

The characteristic function can be given in terms of 
normally ordered moments of the creation and annihilation 
operator ifTsIl : 



V. CONCLUSION 

We have used experimental data sets of quadrature mea- 
surements on phase-diffused squeezed states for a test of dif- 
ferent nonclassicality criteria. Even for a completely phase- 
diffused squeezed state, i.e. where the measured statistics 
were identical for all homodyne detection phase angles, we 
found a pronounced nonclassical character This could be il- 
lustrated with the help of the characteristic function of the P 
function: It directly shows nonclassical features in the lowest- 
order criterion of lfl2ll and the nonclassicality was detected 
with a rather high signal-to-noise ratio. Other nonclassical- 
ity criteria, such as higher-order squeezing or the violation 
of positive semidefiniteness of Agarwal's matrices, fail to re- 
veal nonclassicality beyond squeezing or only show nonclas- 
sical behavior in matrices of higher dimension. Therefore, we 
demonstrated for the radiation under study, that the character- 
istic function of the P function, which contains information 
about all moments of the state, can be a more powerful tool 
for the examination of nonclassicality than a finite set of mo- 
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Introducing the phase-dependent quadrature operator x{ip) = 
e-'^f + ae''^ , we. have 



fc=0 
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(A2) 



Consequently, the normally ordered moments (: x{ip)^ :) can 
be calculated from the characteristic function of the P func- 
tion as 



Qk 



d{ib} 



(A3) 



b=0 



To obtain a formula which can be applied in practice, we insert 
Eq. ( [Tol l into Eq. ( |A3t and use the definition of the Hermite 
polynomials in the form (— l)'^iJfe(^)e^^' = ^''■^ 
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2x2 Matrix 


4x4 Matrix 


6x6 Matrix 


8x8 Matrix 


10 X 10 Matrix 
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6.3 
12.6 
22.2 
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-0.6362(1 ± 0.25%) 
-0.5717(1 ± 0.03%) 
-0.4060(1 ± 0.08%) 
0.0197(1 ± 3.0%) 
1.0000(1 ±0%) 


-4.294(1 ±0.86%) 
-3.337(1 ±0.11%) 

-2.040(1 ± 1.1%) 
-0.2323(1 ± 1.1%) 

0.7856(1 ± 1.2%) 


-104.0(1 ±2.5%) 
-69.93(1 ±0.35%) 

-6.728(1 ±53%) 
-0.5358(1 ± 4.1%) 

0.5493(1 ± 12%) 


-6201(1 ±6.1%) 
-3593(1 ± 0.98%) 
-107.4(1 ± 110%) 
-2.299(1 ± 71%) 
10.85(1 ± 13%) 


-722 ■ 10^(1 ± 12%) 
-335 ■ 10^(1 ± 2.5%) 
-1259 ■ 10^(1 ±49%) 
-459 • 10^(1 ±40%) 
1113(1 ±94%) 



Table 111: Table of minimum eigenvalues of the matrices A/''' for different I and phase noise a. The existence of significantly negative values 
indicate the nonclassicality. 

Hence, we obtain normally ordered moments from measured 
quadratures via 



1 ^ 



The approximation sign indicates that the right hand side is a 
statistical estimator 



^ = z6/ v2. Neglecting the phase argument, we find 
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